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We apply the Implicit Regularization Technique (IR) in a non-abelian gauge theory. We show 
that IR preserves gauge symmetry as encoded in relations between the renormalization constants 
required by the Slavnov- Taylor identities at the one loop level of QCD. Moreover, we show that the 
technique handles divergencies in massive and massless QFT on equal footing. 
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INTRODUCTION AND BASICS OF 
IMPLICIT REGULARIZATION 



Dimensional Regularization (DR) is the natural frame- 
work for computing Fcynman diagrams in gauge field the- 
ories. However the regularization of dimension specific 
quantum field theories such as chiral, topological and su- 
persymmetric gauge theories is known to be a delicate 
matter in the context of DR. That is because the ana- 
lytical continuation on the space-time dimension of the 
Levi-Civita tensor is not well-defined whereas supersym- 
metry is intrinsically defined on the physical dimension 
of the underlying model. 

Although some extensions of DR have been con- 
structed (e.g. Dimensional Reduction pj) they are in 
general inconsistent in arbitrary loop order and may give 
rise to spurious anomalies. Hence a judicious order by 
order calculation in which the symmetry content of the 
model is assured via constraint equations has to be per- 
formed. The drawbacks are clear: in addition to turn- 
ing the calculations cumbersome and tedious, we cannot 
rely on such procedure to study anomalous (quantum 
mechanical) symmetry breaking. This is particularly rel- 
evant in the supersymmetric extensions of the standard 
model 0. 

This motivates the search for a non-dimensional reg- 
ularization/renormalization scheme which, besides pre- 
serving the vital symmetries of the quantum field the- 
oretical model, is friendly from the calculational view- 
point. 

Implicit Regularization (IR) is a momentum space set- 
ting to perform Feynman diagram calculations in regular- 
ization independent fashion. The Lagrangian of the un- 
derlying quantum field theory is not modified: neither an 
explicit regulator is introduced nor the dimensionality of 
the space time is moved away from its physical dimension. 
It has been successfully applied to various quantum field 
theoretical models including those which make sense only 
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in their physical dimension. For quantum eletrodynam- 
ics, theories involving parity violatingobjects (Chern- 
Simons, Chiral Schwinger Model), see [6|. For the study 
of anomalies and CPT violation in an extended chiral 
version of quantum electrodynamics see Q. A compar- 
ison between IR, dimensional regularization, differential 
renormalization and BPHZ forest formula can be found 
in [j| , where it is also calculated the beta function to one 
loop order in quantum chromodynamics. In |Tflj a model 
calculation using (f> 3 theory in 6 dimensions illustrates 
how IR works when overlapping divergencies occur. In 
it is shown that IR is manisfestly supersymmetric in- 
variant. This is illustrated by renormalizing the massless 
Wess-Zumino model and calculating the beta function to 
three loop order. Application to a Gauged Nambu-Jona 
Lasinio model can be found in 

The main idea behind IR is very simple. The ultravi- 
olet behaviour of the amplitude is isolated as irreducible 
loop integrals (ILI's) which are independent of the ex- 
ternal momenta and need not be explicitly evaluated to 
display the physical content of such amplitude. This can 
be achieved by judiciously using the identity at the level 
of the integrand 
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in order to eliminate the external momentum ki from the 
ILI, N being chosen so that the last term is finite under 
integration over k. 

You may assume, to be very strict, that a regulariza- 
tion (say dimensional regularization) implicitly acts on 
the amplitude in order to use Q in the integrand. How- 
ever once you have separated the divergencies as irre- 
ducible loop integrals from the finite part of the ampli- 
tude you need not compute the divergent integrals within 
IR. They may be subtracted and absorbed in the coun- 
terterms exactly as they stand. The explicit computation 
of such ILI's is the origin of spurious symmetry break- 
ing which may contaminate the physics of the underlying 
model. In [g we have defined what is meant by a min- 
imal subtraction scheme within IR and compared with 
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dimensional regularization and differential renormaliza- 
tion. In this process a natural renormalization group 
scale emerges as it should. The generalization of this pro- 
gram to higher loop order is straightforward: the over- 
lapping divergences can be treated in a similar fashion 
following a well-defined pre scrip tion which corresponds 
to the BPHZ forest formula [ljf . 

At the one loop level in the Minkowskian 
4— dimensional space time the ILFs show up as 
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where J fc = J (f^yr , k being the internal momentum, /i 
is an infrared regulator and n = 1,2, .... Typical higher 
loops (logarithmically divergent) ILI's are, in 4 dimen- 
sions 
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where 4 + m = 2p and A is an arbitrary dimensionful 
non-vanishing constant originated in the previous order 
0. Some comments are in order. For massless models 
we may always introduce a ficticious mass to regulate 
the propagators in the infrared limit without sacrificing 
neither gauge symmetry Q nor supersymmetry Q. We 
shall explicitly verify this in the context of QCD. Local 
arbitrary counterterms will appear in IR as differences 
between irreducible loop integrals of the same degree of 
divergence. Because we are not explicitly evaluating the 
divergent integrals, such (finite) differences will have the 
status of free parameters which should be adjusted by 
phenomenology or symmetry constraints. Explicit regu- 
larizations will generally assign a (regularization depen- 
dent) value to such differences which may lead to sym- 
metry breaking. 

In 3 + 1 dimensions at the one loop level these arbitrary 
parameters look like: 
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g{^g a p} stands for g^g a p + g m g v p + g^g va , and the 
ctj's are arbitrary, finite and regularization dependent. 
Similar relations appear at higher loop order. 

It is straightforward to see that in dimensional regu- 
larization (0J-® evaluate to zero. In |(| we have shown 
that vector gauge symmetry is compatible with setting 
all the ai's to vanish. However we have shown that this 
is not the only solution. Such feature somewhat explains 
why dimensional regularization is gauge invariant. 

Whereas fixing the a,'s to zero from the start is more 
practical from the calculational viewpoint care must be 
exercised when dealing with dimension specific objects 
such as axial vertices and Levi-Civita tensors, such arbi- 
trary parameters should be fixed on physical grounds. 
In [6j we demonstrated that are connected to 

momentum routing invariance in a Feynman diagram. 
Should on vanish then the amplitude is momentum rout- 
ing invariant. The ideal arena to test this feature is 
the study of chiral anomalies: in perturbation theory 
such anomalies manifest themselves as a breaking of mo- 
mentum routing invariance In we have stud- 
ied the Adler-Bardeen-Bell-Jackiw anomaly for arbitrary 
momentum routing and seen that IR consistently display 
the triangle chiral anomaly in a scheme-free fashion in 
a way that the anomaly appears in the vector and ax- 
ial Ward identities on equal footing. This is the best 
we can expect from a regularization scheme Q . We have 
also seen that the mass spectrum of the Chiral Schwinger 
model is undetermined by a arbitrary paremeter which 
in perturbation theory corresponds to a finite arbitrary 
number expressed by the difference between logarithmi- 
cally divergent integrals in IR. This is what is expected 
from a non-perturbative calculation Q. 

In other words we have seen that should such differ- 
ences be set to zero {04 ' s — 0) then the amplitude is mo- 
mentum routing invariant and (abelian) gauge invariant 
(although if they assume a non-vanishing value it does 
not necessarily mean that gauge invariance is broken). 
When an explicit form of the regulator is used such dif- 
ferences are assigned a (regularization dependent) value. 
In general we should keep any arbitrariness which ap- 
pears in perturbation theory until the final stage of the 
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calculations where physical conditions may fix its value. 
In this sense, IR is especially taylored to implement this 
idea especially when quantum symmetry breakings may 
occur. 

In Q we have verified that constraining to zero also 
ensures the transversality of the vacuum polarization ten- 
sor of QCD. The next stringent test for establishing the 
generality of IR is to extend these ideas to nonabelian 
gauge field theories. 

The purpose of this paper is threefold: 

1) To check whether a constrained version of IR (CIR) 
generalizes to a nonabelian gauge theory (QCD) and 
show that gauge symmetry is preserved as expressed 
by the Slavnov- Taylor identities between the renormal- 
ization constants and calculate all the renormalization 
group constants to one loop order; 

2) To define a renormalization group scale in the fol- 
lowing way: just as in the case of the photon propagator, 
we introduce a ficticious mass (p.) for the gluon which will 
appear in both the finite and the (logarithmically) diver- 
gent pieces of the amplitude. For the divergent piece it 
will show up as the ILI 
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We eliminate the infrared mass regulator from the defi- 
nition of the counterterm by using the identity (see ap- 
pendix A) 
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with D ab and D l J referring to the adjoint and fundamen- 
tal representation of the colour group SU(3). Also 



a is the gauge-fixing parameter and A 0fi are the gauge 
fields coupled to n/ Dirac fermions ipo and to the ghost 
fields Co- The index "0" stands for bare quantities. The 
group theoretical factors which will appear in the am- 
plitudes are defined through the relations tr(t^t b ) = 
C(r)S ab , t a r t a r = C 2 (r)i, f acd f bcd = C 2 (G)S ab . Because 
the interaction terms in the Lagrangian above are inter- 
related by BRS symmetry, only one coupling constant is 
left independent. Consequently the renormalization con- 
stants will be constrained by generalized Ward-Takahashi 
identities ( Slavnov- Taylor identities) . 

We define renormalized fields and couplings through 
the renormalization constants as below 
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and A is a nonvanishing parameter which parametrizes 
the arbitrariness in separating the divergent from the fi- 
nite content of the amplitude and will play the role of 
a renormalization group scale in IR. As a byproduct we 
realize what is meant by a minimal subtraction, mass in- 
dependent scheme in IR, namely subtracting Ii og {\ 2 ) . 
The infrared divergent term expressed by In /i 2 will ex- 
actly cancel the infrared cutoff dependence in the finite 
part of the amplitude, as it should, for all infrared safe 
theories. 

3) to see that unlike Dimensional Regularization, the 
tadpole graphs of Yang-Mills fields play a crucial role for 
maintaining manifest gauge invariance through cancella- 
tions of quadratic divergences which appear at one loop 
order. 
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Therefore we may define Cq = C + C c t where C is pre- 
cisely equal to Cq except that it is written in terms of 
the renormalized variables whereas C c t is the countert- 
erm Lagrangian which reads 

Cct = {Zz-l) l -A»5 a \g^d 2 -d^d u )A 1 (,+ 
+ (Z 3 - l)c a 5 ab (-id 2 )c h 

- (Z 2 Z m -l)m^ 1 

- {Z 1 - V-gf^id^ -d„A*)AZA» 

(z 4 -i)i 9 2 r b 7 cde A Q X^ 

- {Z t - l)igf abc (d»c?)c b A^ 
+ (Z 1F - 1)9$%^ 'A% , 

where we have defined 
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The QCD Lagrangian reads 
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The equality of Z g for all the couplings leads to the 
Slavnov- Taylor identities: 



Z 3 



Zi 



Zip 
Z 2 



Zi 
Zx 



(15) 



The Feynman rules for QCD can be found in any text- 
book. We follow T. Muta and work in the Feynman 
gauge where a = 1. To one loop order, the relevant am- 
plitudes are represented by well-known diagrams which 
we shall call as 
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V ab 
ghost 

K 

a abc 

a abc 

P- 

a abed 



gluon self-energy 
quark self-energy 
ghost self-energy 
quark-gluon vertex 
three-gluon vertex 
ghost-gluon vertex 
four-gluon vertex . 



(16) 



We start with the gluon self-energy which is composed 
of four contributions as depicted in fig. Q). 

K = n #(!) + n ^( 2 ) + K(V + ( 17 ) 

where n°»(l), IIJ* (2), IIJJ*(3) and nj*(4) represent the 
quark loop, the gluon loop, the gluon tadpole and the 
ghost loop respectively. It is purely transversal as re- 
quired by the Slavnov- Taylor identities and thus it does 
not admit a mass term and there should be no mass 
renormalization. Hence the quadratic divergences which 
appear in 11"^ should cancel out. 



(1 1 01) may be cast as 

n#(2) = -\g 2 C 2 {G)5 ab [{2p p p„-Ap 2 g pu )J(p 2 ^ 2 ) 
- g^(2I q uad(p 2 )+P a P^%) 
- W( Pl/ J^p 2 ^ 2 ) + J^ip 2 ^ 2 ))}. (22) 



As for the ghost loop, we have 



n£(3) = -g 2 f dac f cbd 



= -g 2 S ab C 2 (G)( Pu J p (p 2 ^ 2 
+ J^{p 2 ^ 2 )) , 



in which J M „, J a and J are defined as 
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(25) 
(26) 



Mp 2 ,H 2 ) = -2p Q e(°) + -p,Z{p 2 ^ 2 ), (27) 



J = I log ( f i 2 )-bZ(p 2 , f i 2 ), 



with 
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which, for — > 0, is given by 



IC(1) = -g z C 2 {G)5 ab i 



9^v 



k k2 ~ M 2 

- -ig 2 g^C 2 {G)5 ab I quad {^ 2 ) (18) 

where p, is an infrared cutoff (mass regulator) which 
should be set to zero in the end. 
The gluon loop amplitude reads 



U ab (2) — - I 2 f acd f bcd M . 



where 
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N tw _ \gi>P(p-k) a +g pa {2k+pY +g <T ' M (-k-2p)P] 

x [S; (k-p)„ + g pa {-2k - P y + 6»(k + 2p) p ] 

= 2p ll p v - 5{p^,k v + p v k p ) - lOk^K 

- g^[(p~kf + {k + 2p) 2 ]. (20) 

Using that 

(p - k f + (k + 2p) 2 = {k+pf + k 2 +4p 2 , (21) 



In 



e 2 [i 2 



(30) 



Collecting all the results so far enables us to write 
£n£(i) = g 2 C 2 (G)S ab \ (pV~ W) x 



+T 2 +» 2 T (0) +» Q »' 3 T (0) + 



+P + p^T^ + p , (31) 



in which the T's are the arbitrary constants defined in 
the relations (0J to Q . Moreover in writing l|31|l we have 
absorbed some constant factors in the T's. 

The fermion loop contribution to the gluon self energy 
is identical to the vacuum polarization tensor of QED ex- 
cept for the colour and number of fermions (n/) factors. 



5 



It has been computed within IR elsewhere 0. With- 
out loss of generality we write the result in the limit of 
massless fermions to yield 

D#(4) - p 2 C(r)n f S ab ^p,p„ -p 2 g^) x 



, T 2 , D 2 T (0) , D a„/3 T (0) , 



- fXV 



+PV$1 +pPpuT { $ +pV^T$ j • (32) 

Some comments are in order. Firstly notice that the 
quadratic divergences expressed by g^Iq U ad{^ 2 ) and 2 
that appear in the gluon tadpole, the gluon loop and the 
ghost loop amplitudes combine to make up Y 2 „ = a\g^ v . 
Gauge invariance tells us that we ought to set ot\ — as 
well as all the others CKj's as defined in JH-Q. Hence 
tadpole graphs of gauge fields play an essential role in 
maintaining gauge invariance within our framework. DR 
automatically sets quadratic divergences to zero in the 
limit where fj, — ► 0. Here this is not necessary in order 
to ensure the transverse form of the gluon self-energy as 
required by gauge invariance. As we shall see, setting 
Aj's to zero in Q to (7J automatically preserves (vector) 
gauge invariance through the Slavnov- Taylor identities. 
This is in consonance with the idea that ultimately one 
should let arbitrary parameters to be fixed on physical 
grounds. In the present case, gauge invariance does this 
job. However they were shown to play a crucial role 
in describing correctly chiral field theories in which the 
Dirac algebra involving 7^ matrices prevents the use of 
naive DR. In recent work [12j , such free parameters have 
been taken into account in a renormalized version of a 
SU(3) Nambu and Jona-Lasinio like model. The rele- 
vant observables have been calculated in excelent agree- 
ment with experiment including a simultaneous and sat- 
isfactory fit for both /„■ and f K . This is an interesting 
feature because it enables IR to be applicable to study 
the dynamics of effective field theories (for instance the 
derivation of the gap equation in the gauged Nambu and 
Jona-Lasinio model [13j in a gauge invariant fashion). 
Moreover the leading quadratic terms also play a crucial 
role in the evaluation of the hadronic matrix elements of 
four quark operators in the kaon decays K — ■> 7T7t as well 
as in providing a consistent prediction on the direct CP 
violating parameter e'/e in kaon decays |l4j . IR may be 
applied to all these scenarios. Operationally it is conve- 
nient because one has a gauge invariant momentum space 
framework. 

Note that the algebraic procedure that we have used 
to define a mass independent, minimal, renormalization 
scheme naturally introduces an arbitrary scale A. As we 
shall see A plays the role of a renormalization group scale. 

In order to define genuine renormalization constants 
which display the ultraviolet scaling behaviour of the 



model we use the identity © and note that the (infrared) 
divergence parametrized by In /i 2 as /1 — > cancels out 
against an identical term coming from the UV finite part 
whilst an arbitrary nonvanishing parameter A appears. 
Altogether m* = £i=i nj* (i) reads 
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A 2 



+ (15C 2 (r) -6n/)ln(-j 



2C 2 (r) + 2n/j 
(2 3 -l)^(p 2 V-wJ- (33) 



We define the counterterm for the amplitude (|33() by 
minimally subtracting (in the IR sense) the ILI expressed 
by Iiog(^ 2 ) to define 



; C 2 (G) - -n f C(r) 



/ log (A V + 0( 5 3 ). 



(34) 

Note that the algebraic procedure that we have used 
to define a mass independent, minimal, renormalization 
scheme naturally introduces an arbitrary scale A. As we 
shall see A plays the role of a renormalization group scale. 
In [8| we compare renormalization schemes in IR, DR and 
differential renormalization (see also (lif). 

The quark self energy T,(p) is similar to the electron 
self energy apart from a group theoretical factor. £(p) 
has been calculated in Q within IR. It reads 

E(p) = -ig 2 C 2 {r)(jp- 4m)/ log (m 2 ) + g 2 g^T^p 
+ (Z 2 -l)p-{Z 2 Z m -l)m + ^(p,m), (35) 

where the tilde means that the quantity is finite. We 
shall use such notation from now on. Indeed S(p, m) is 
both ultraviolet and infrared finite. In order to define the 
corresponding renormalization constants we consistently 
make use of (0 in order to pursuit a mass independent 
scheme as well as introducing the arbitrary constant A 2 . 
It is noteworthy that the b In fj, 2 piece coming from © 
cancels exactly the infrared divergence in the ultraviolet 
finite portion of the amplitude, as it should. 

Henceforth we shall systematically define the renor- 
malization constants in a mass independent fashion 
(which defines the "minimal" scheme in IR) as well as 
set a-i's (T's) to zero. We define this procedure as con- 
strained IR (CIR). We shall no longer write the T's 
explicitly in the remaining amplitudes for the sake of 
brevity. 

Therefore the fermion mass and field renormalization 
constants can be cast as 



Z m = l + 3zC 2 (r)/ log (AV + 0(g 3 ) 
Z 2 = l + iC 2 (r)I log (\ 2 )g 2 + 0(g 3 ). 



(36) 
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FIG. 1: 1 loop diagrams of QCD 



After the appropriate color index contractions, the one 
loop correction for the ghost propagator simplifies to 

KhostiP 2 ) = 9 2 9^C 2 (G)5 ab ^ [{p - ^ 2 k J 2] P {k2 

+ (Z 3 -l)S ab P 2 , (37) 

where fi 2 is an infrared mass regulator for both the ghost 
and gluon propagators. After some straightforward alge- 
bra we have 

KlstiP 2 ^ 2 ) = ^V(^C 2 (G)/ log (A 2 ) 

+ Z 3 -l)+Z a g b host (p 2 ,\ 2 ), (38) 

from which we define the renormalization constant 

Z 3 = l- l -C 2 {G)I Xog (\ 2 )g 2 + 0(g 3 ). (39) 

For the one loop quark-gluon vertex A^ shown in fig. 
(J2J) we have two contributions: the QED-like electron- 
photon vertex diagram A®(p, q)(l) and the one involving 




FIG. 2: 1 loop diagrams of QCD 

the three-gluon vertex A"(p, q)(2). The former differs 
from the QED electron-photon vertex by a group theo- 
retical factor 

AJ(p J? )(l) = (iW)AjP D (p,«). (40) 

We have also computed A^ ED (p, q) within IR in [8j so 
here we only quote the result: 

iA% ED (p, q) = 7 >2 + Wm 2 )] + A"(p, q) , (41) 

where m is the mass of the^Termion, A M (p, q) is finite and 
a 2 is arbitrary and shall be set to zero within constrained 
IR. Using that t d t a t d = [C 2 {r) - 1/2 C 2 (G)]t a enables us 
to write 

A>, g ,A 2 )(l) = -ig 3 t a ^(c 2 (r) - \c 2 {G)) x 
x(l log (\ 2 ) + A a (p,q,\ 2 )(l)). (42) 
As for A"(p, q) (2), the Feynman rules give 

l A«(p 7 q)(2)=g 3 r bc t b t c Jj± (43) 
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with 

K = 'f(/t+myf((2p-q-k) v g w> 

+ (2k-p-q) p g fMU + (2q-p-k) fl g U pJ , 
V = (k 2 - m 2 )[(k - pf - p 2 ][(k - qf - p 2 ] . (44) 

where again p is a mass regulator for the gluon propaga- 
tor. We proceed as before. We remove the external mo- 
mentum dependence from the ILI by applying the iden- 
tity Q in the propagators which contain the momenta p 
and q above. Then we isolate a genuine (ultraviolet di- 
vergent only) contribution for the counterterm with the 
help of identity |J5J| to get, in CIR, 

A>, g ,A 2 )(2) = ^g 3 C 2 (G)t a ^I log (\ 2 ) 



+ A>,g,A 2 )(2) 



(45) 



Finally we define the renormalization constant Z\p to 
one loop order by adding up the two contributions for 
the quark-gluon vertex: 

A>,<?,A 2 ) = g~ ftl t a (-ig 2 (C 2 (G) + C 2 (r))I log (\ 2 ) + 
+ Z 1F -l + A a {pq,X 2 )^ (46) 

which gives 

Z 1F = 1 + ig 2 (c 2 {G) + C 2 (r))l log (X 2 ) . (47) 

To calculate the renormalization constant Z\ we work 
with the ghost-gluon vertex. It receives two contributions 
(fig. |2J): the ghost-ghost-gluon loop A^ c and the ghost- 
gluon-gluon loop A^f. Let pi, p 2 and p\ +p 2 = q be the 
external momenta. Thus 

ibc X a-bc i a afc. 



K c = K c + K c + H)gf abc ^(Zi - 1) 



(48) 



It is straightforward to use the Feynman rules with the 
help of the identity f-fefepfcpf = N/2f abc for SU(N) 
to get 



A^ c (Pi,P2) = -yC 2 (G)fV 



k (k2-p2)[(k- Pl )2-p2][(k-qr-p2] 



(49) 



We proceed according to the rules of CIR, as we have 
done before, to arrive at 



A?>i, P2 ) = ^C 2 (G)r bc q,I log (\ 2 )+A?;(p llP2 ) 



Similarly we have for the other contribution 



■iff 



(50) 



A^ c bi,P2) = -j-C 2 (G)r bc qfl I log (\ 2 ) + A£ c (p 1)P2 ) 



and hence 

kf c {puP2) = - J .gr hc g Al (-z^C 2 (G)/ io9 (A 2 ) + Z 1 -l) 
+ ~Af c (pi, P2 ). (52) 

Finally wc define the renormalization constant in a min- 
imal fashion within IR as 



Zi 



l + ^g 2 C 2 (G)I log (X 2 ). 



(53) 



The class of one loop three-gluon vertex graphs from 
which we shall define Z\ are shown in fig. ^ They have 
been explicitly computed in [l^, 01 within DR. The 
calculation is straightforward yet tedious. We have pro- 
ceeded according to the rules of CIR as before in order to 
isolate the ultraviolet divergence as a term proportional 
to Ii og (X 2 ) after making use of JJjJ. The infrared diver- 
gent piece proportional to ln(A 2 /^ 2 ) as p — » cancels out 
with an alike term stemming from the ultraviolet finite 
piece of the amplitude as we generically discuss in the 
appendix B. For the sake of brevity we shall present only 
the result here. Let p and q be the external momenta. 
Then 

Kl c x (P, <z) = ~*gf abc V^x (p, q, p + q) x 



-ig^-lc 2 (G) + ^C(r)n f )l log (X 2 ) 



+ Z l -l\+KH\{pq) 



(54) 



(51) 



V^x{p,q,p+q) = {p-q)xg )J ,u-p^guX+qug^x, from which 
we define 

Zi = 1 + ig 2 ( - \c 2 {G) + ^C(r)n f j I log {\ 2 ) . (55) 

Last but not least we have to compute the four gluon 
vertices depicted in fig. [5] along with all their permu- 
tations. This long calculation has been performed with 
great detail by Pascual and Tarrach in [TtI in the Wein- 
berg's scheme as well as by Papavassiliou in using 
the S-matrix pinch technique . The corresponding result 
within CIR reads: 

a$£,(pi,p2,ps,P4) = -g 2 wXu x 

X ( ~ % ( C ' (G) + AC ^ n f) h °3 ( A ') + Zi - 

+ %$L(pi,P2,Pa,Pi), (56) 

where we have used the same notation as T. Muta [Tlj. 
namely 

T / t/a 1 a2a3a 4 _ f j- 13,24 _ j-14,32\ 

+ (/ 12 ' 34 - Z 14 ' 2 J.9MiM3.9m2M4 

+ (/ 13 ' 42 -/ 12 ' 34 ).9, 1W 5M3, 2 , (57) 
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with pi' 1 ™ = f^ ] ajaa l a m _ Therefore 



Now using that 



Z 4 = l + ig 2 f hj 2 {G) + |c(r)n, 1 h og {X 2 ) . (58) 



III. SLAVNOV- TAYLOR IDENTITIES AND 
RENORMALIZATION GROUP FUNCTIONS 

It is a simple task to verify that CIR explicitly pre- 
serves the Slavnov- Taylor identities expressed by (|15fl : 

| i = | 1 = ^ = | i = l + ^ 2 (G)/ ;o3 (A 2 ). (59) 



Z. 



In other words CIR explicitly fixes the arbitrariness of 
IR in such a way that gauge invariance is maintained. 
At higher loop order similar relations to those displayed 
in equations Q-Q are expected to hold Q and its con- 
strained version should implement vector gauge invari- 
ance as well poj ]. 

In defining a mass-independent minimal scheme in CIR 
there appeared an arbitrary non- vanishing constant A 2 . 
As discussed before, by subtracting only the term propor- 
tional to Ii og (X 2 ) defines a minimal subtraction scheme 
within IR and we are left with the finite piece of the 
amplitude which is also dependent upon A 2 . Moreover 
it is identical to the amplitude that we would obtain 
had we employed differential renormalization [T^ |. 
The arbitrary scales which appear in IR and Differential 
renormalization can be identified and thus the truncated 
connected n-point renormalized Green's function, say 
G^\pi, g,m) is expected to satisfy a Callan-Symanzik 
like renormalization group equation where A plays the 
role of renormalization group scale. Thus we have 

( d d d 



-n A jA(g)-n m (g) G« =0, 



(60) 



where ua {nf) is the number of gluon (quark) legs in 
momentum space, m and g are defined as in 1131) . 7a(7i/>) 
are the anomalous dimension of the gluon (quark) field 
and 

m = a* 



im{g) 
iA{g) 

7V>(sO 



A dm 

m dX 
A dZ 3 

2Zid\ 
X dZ 2 



and 



(61) 



2Z 2 dX ' 

For instance let us explicitly calculate the (3 function. 



Recall (JT3J): g = Z g g, Z g = Z ± Z 3 
d 



3/2 



Hence 



2X 1 



dX- 



\z g g) 



(62) 



(63) 



in the equation above yields after some simple algebra 

3 

/8 = -3^3a(llC 2 (G)-4C(r)n / ) +0(g 5 ). (64) 

In a similar fashion we may use the renormalization con- 
stants which we have calculated in CIR to show that 



6g 2 
(4tt) : 



;C 2 (r)+0(g 5 ) 



(65) 



1A 



3(4tt) : 



-(bC 2 {G) - AC{r) nf ) + Q(g 5 ) , (66) 



and 



7V 



(47T) 



;C 2 (r)+0(g 5 ) 



(67) 



which are the standard values of the renormalization 
group functions. Particularly in a minimal scheme within 
our framework, they coincide with the MS scheme in di- 
mensional regularization. 



IV. CONCLUSIONS 

We have applied the IR method to QCD at the one 
loop level. We have shown that it preserves non abelian 
gauge symmetry and that there is no need of a different 
prescription to deal with massive or massless theories. 
A constrained version of IR which implies in momentum 
routing invariance also delivers gauge invariant ampli- 
tudes for the non-abelian case. 



Appendix A: 

Consider a bubble (or a piece of a certain QCD ampli- 
tude) : 



I" 



A 



4-d 



1 



(fc 2 — [i 2 )[(k + p) 2 — [i 2 



(68) 



The index d stands for dimensional regularization. Al- 
though IR docs not use a explicit regulator, we will use di- 
mensional regularization here for pedagogical purposes to 
show how to define a irreducible loop integral in a mass- 
less theory free of infrared divergencies through equation 
(EJ . Similar relations can be derived at higher loop order 
I 

We follow the prescription of IR to separate the diver- 
gent parts by means of the identity given by equation 
(|T|) . As we have discussed we may strictly assume an im- 
plicit regulator to manipulate algebraically the integrand. 
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However as we do not actually evaluate the irreducible 
loop integrals we need not to make a regulator explicit. 



1 



k (fc 2 - M 2 ) 2 
p 2 + 2p ■ k 



k (fc2- M 2 )2[(A . +p)2 _ M 2] 



(69) 



with e = 2 — d/2. As a matter of illustration we calculate 
the first integral (Ii og (m 2 )) using dimensional regulariza- 
tion to obtain 



k 2 'I log (m 2 )=b 



- + A + In 

e 



4A 2 
m 2 



+ 0(e), (70) 



where A is a constant characteristic of dimensional regu- 
larization. The second integral is finite and evaluates to 



rd , f\ , ^ 2 Z(1-Z)-/X 2 

V = -^ o dzlnj^ — 2 

In the limit where /i 2 — > 0, we have 



rd 



-6 In 



e 2 /i 2 



In IR we write 



/ d =I to5 ( M 2 )-Mn(- 
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However because 
A 2e / io9 (A 2 )-6 
A / we have 



A + ln 



4A 2 
A 2 



A^ os (^)-A^J Joff (A^) = &ln l , 

v A* 



(71) 



(72) 



(73) 



+ 0(e), (74) 



(75) 



which is just equation @ in the limit e — > 0. Substituting 
this relation in the expression for I d yields 



7 d = / /o9 (A 2 )-61n 



e 2 A : 



(76) 



Now we are allowed to subtract a genuine ultraviolet di- 
vergent object Ii og {\ 2 ) by defining the appropriate coun- 
tcrtcrm. The non-vanishing arbitrary parameter A plays 
the role of renormalization group scale. 



Appendix B: Infrared Finiteness of the One-Loop 
Amplitudes 

We now turn ourselves to an important discussion on a 
problem which arises when the limit p? — > is taken. We 
must make sure that, when using the scale relation given 



by equation (7), the term in ln(/x 2 /A 2 ) will be cancelled 
out by a contribution that comes from the ultraviolet 
finite part. The divergent integrals that are present in the 
calculations at the one loop level for the renormalization 
of QCD are 



.4 



1 



.4 
B 



k (k 2 -^ 2 )[(k + p) 2 -p 2 y 

kit 



" ~ ' k (k 2 -^ 2 )[{k + p) 2 -p 2 Y 

» V ~ 'k (k 2 -^ 2 )[(k+p) 2 -p 2 }[(k + q) 2 -p 2 \ 



(77) 
(78) 
79) 



k (k 2 - v 2 )[(k+p) 2 - p 2 ][(k + q) 2 - p 2 \ 



■80) 



In the integrals above we have introduced the mass /j, 
that will be set to zero in the end. Since the integrals are 
assumed to be regularized, we can follow the prescrip- 
tion of Implicit Regularization, and separate the diver- 
gent parts by means of the identity given by equation (1). 
Below we show the calculations: 



• A: 



A = 



1 



k (fc 2 -M 2 ) 2 

p 2 + 2p ■ k 
k (k 2 - fi 2 ) 2 [(k+p) 2 - fi 2 



(81) 



We now use the identity expressed by equation (7) 
in the first integral. The second one, which is finite 
and does not depend on any specific technique, is 
given by 



, f 1 , . fp 2 z(l-z)-fi 2 



For n 2 — ► 0, we have 



A fi n = b In ( - 



2 2 



(82) 



(83) 



We clearly see that the /i 2 cancels out when the two 
parts are put together, so that we obtain 

-2 



A = I log (\ 2 ) - bin 



P 



3 2 A 2 J ■ 



(84) 



• A^: After the expansion, some integrals vanish, 
since their integrands are odd in the integration 
variable. After calculating the finite part, we have 



An = -2p v 



(k 2 - /x 2 ) 3 



b ( p 



e 2 fi 2 



(85) 



The divergent integral is Q$([A 2 ) (see eq. (6)) and 
we use equation (3) to write 

eftV) = 9 -f (Wm 2 ) - a 2 ) 



^^I log (X 2 ) + b\n^-)-a 2 ) (86) 
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Again we see the cancellation of fj, 2 when the finite 
part is considered. We are left with 



.4 



Pu- 



e 2 \< 



Ct2 



(87) 



It is important to note that equation (3) is essential 
in the cancelation of the y 2 . 

• B^ v : After the expansion and calculation of the 
finite part, we obtain 



Buv — 



1 



In 



P 



e 2 \j? 



2 \ Q 2 P 2 
2 + M Voo 1 - y'ho - y Voi 

p IJ .PuVo2 + q^qvmo + ip^qv + q^Pu)mi} >' 



where 



dz 

I! JO 



dy 



z n y m 



Q(p, q,y,z) 



(89) 



and 



Q(p,q,y,z) = p 2 y(i-y) + q 2 z(i-z) 

- 2(p- q)yz - 



(90) 



It can be easily seen that the functions rj nm do not 
have problems when /i 2 — » 0. The divergent inte- 
gral has the result of equation (|86[1 , so that 



B„ 



9 -f[l log {\ 2 ) bin 



P 



Oil 



+ f(Vnm), (91) 

where f(r) nm ) represents the r\ dependent part. 
• B^ va : The mechanism is the same as in the other 



integrals: 

B, va - -2(p + q) J k{k2 _ fl2)4 

+ To" [(P + Q)n9ua + {P + q)v9va 



+ (p + q) a g^ v ] in , ., 



+ 9(Vnm)- (92) 



The 5(77) represents the r\ dependent part. The 
divergent integral is ©^^(a* 2 )? a s defined in equa- 
tion (6) . We use the relation given by equation (5) 
to write 



1 



24 

+ felll ( ) n j 



(93) 



The substitution of this result in B^ va and the 
adoption of the same procedures as before, lead us 
to the cancellation of the infrared divergences and 
we have the final expression 

B^ a = -—[(p + q)ng va + (p + q)vg m 

+ (p + q) ] (WA 2 ) 



fehl - 



2 ? 

e [i 



- a 4 ) + g(r] nm ). (94) 



We call the reader's attention to the fact that, in the last 
calculation, the equation (5) was mandatory in order to 
/1 2 to be cancelled. It is also interesting to note that 
the same relations that are necessary to preserve gauge 
invariance are also essential for the cancellation of the 
infrared cut-off /i 2 . 
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